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mental standards, and by statistical and systematic 
errors in the measurement of x-ray lines and in the 
separation of line intensity from spectral back- 
ground. The errors in estimating the characteristics 
of the detector system, with exception of deadtime 
effects, cancel when the emission from the speci- 
men is divided by that from the standard. The esti- 
mate of achievable accuracy of EPMA and the 
comparison of diverse approaches must therefore 
involve consideration of theoretical aspects, of er- 
rors in the values of parameters used in the proce- 
dure, and of operational errors in specimen and 
standard preparation and in analysis. In the devel- 
opment of current "correction procedures," theo- 
retical, mathematical and measurements aspects are 
interwoven. For instance, the distribution in depth 
of x-ray generation can either be calculated by 
Monte Carlo calculations, which are limited by the 
accuracy of available parameters, or by exeriments 
with sandwich tracer targets. To further compli- 
cate things, some Monte Carlo calculations use 
techniques empirically adjusted to fit the available 
experimental evidence. Since the precision of mea- 
surement of relative x-ray intensities is better than 
the accuracy of results on specimens of known 
composition, there is room for improvement; but, 
to achieve this, the diverse potential source of er- 
rors will have to be unraveled and tested sepa- 
rately. 
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1. Introduction 

Secondary Ion Mass Spectrometry [SIMS] pro- 
vides a unique analytical potential due to the fact 
that all elements can be analyzed with high sensiti- 
vity and isotopic specificity, that elemental and 
molecular (bonding) information can be gained and 
that micro, surface and bulk analysis can be perfor- 
med. The major limitations of SIMS are deter- 
mined by the complexity of the mass spectra, the 
large variation of the secondary ion yields for dif- 
ferent elements, and also for a particular element in 
different matrices (chemical matrix effect), and by 
the disorder induced into the analytical zone due to 
the high energy (E =l— 20keV) of the primary 
ions [1,2]. 

To overcome these limitations in praxi, the fol- 
lowing approach for quantitative trace analysis is 
most useful: 

• Use of reactive ion sputtering (oxygen or ce- 
sium) to increase the secondary ion yield and to 
achieve a chemical modification of the surface 
zone which reduces the chemical matrix effects. 

• Use of high performance instrumentation 
which allows the identification (and separation) of 
interfering species by high mass resolution. 

• Quantitation with relative sensitivity factors 
(RSFs) obtained from external or internal calibra- 
tion. 

• Systematic study of all sources of analytical 
errors, particularly of the measurement process, as 
described by Powell [3], 

• Development of problem oriented measure- 
ment techniques and analytical strategies, combin- 
ing SIMS with other analytical techniques yielding 
confirmatory and supplementary information (ex- 
ample ref. [4]). 

Quantitation procedures and the general prob- 
lems associated with quantitation of secondary ion 
mass spectra are discussed in several excellent pub- 
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lications [1,2,5-7]; thus this paper will concentrate 
on specific tasks of quantitative analysis with 
SIMS. These tasks are: multielement, ultratrace 
analysis of metals; high-accuracy, depth distribu- 
tion analysis of dopants in silicon; and quantitative 
interface analysis. 

2. Multielement Ultratrace Analysis of 
Metals 

The need arises mainly from microelectronics, 
with its extreme purity requirements for metals 
used for the production of electrode lines and inter- 
connects in VLSI devices: In refractory metals 
used for 1 Mbit and 4 Mbit storage devices, alkali 
elements and U, Th must not exceed 10 ng/g; the 
concentration of transition metals must be below 
100 ng/g. 

Due to the limited potential of NAA for ultra- 
trace analysis of refractory metals, only mass spec- 
trometric techniques can be considered for direct 
solid state analysis. A comparison of the useful 
yields [8] (SSMS : 10- 8 -10 7 , GDMS 10" 7 , SNMS 
10" 9 , SIMS lO-'-lO- 5 [using reactive ion bom- 
bardment]) shows the large basic potential of 
SIMS. 

The need to establish RSFs with reference mate- 
rials characterized for trace elements — a disad- 
vantage compared to SSMS, GDMS and SNMS— 
makes the analytical procedure more elaborate. 
Figure 1 contains the basic steps for the prepara- 
tion and characterization of the reference material 
(doping with 40 elements, in the ug/g-range, round 
robin analysis, test of homogeneity with SIMS) and 
the analysis of the materials [9]. 

Thirty elements are analyzed with oxygen pri- 
mary ions, 10 with cesium. Interferences are elimi- 
nated by energy filtering, peak stripping and high 
mass resolution. Several mass spectra or depth pro- 
files are recorded to obtain a fairly representative 
value within a reasonable time (ca. 2 hours). Mate- 
rial consumption per analysis is ca. 1 jig (1-10 mg 
for SSMS, 1-10 mg for GDMS). With optimized 
measurement techniques, detection limits (calcu- 
lated from the analysis of the reference materials) 
between 1 pg/g and 100 ng/g are obtained with 
SIMS (SSMS 10-100 ng/g, GDMS 1-10 ng/g) 
[9-11]. 

For the assessment of accuracy, chemical matrix 
effects were studied by comparing results obtained 
with oxygen bombardment (with and without satu- 
ration) and cesium bombardment and using differ- 
ent external standards (doped materials and ion 



implants). The values lie typically within a factor 
of 2 [10]. Furthermore, a comparison of SIMS with 
GDMS and SSMS (including NAA results on a 
few elements) for W and Mo showed that the 
agreement of SIMS with the other methods is on 
the same order (typically within a factor 2) as be- 
tween GDMS, SSMS and NAA, respectively. This 
indicates the absence of large systematic errors in 
the SIMS results due to a chemical matrix effect. It 
was found that the major source for deviations be- 
tween the different methods is the inhomogeneity 
of the materials analyzed. This again shows the 
need for reference materials that are homogeneous 
on the nm to cm scale. 

For homogeneous materials, an analytical accu- 
racy of SIMS for the ultratrace analysis of U and 
Th in aluminum of 10% has been reported [12], 
This corresponds well with the experiences ob- 
tained in the analysis of semiconductors. 

With these figures of merit, SIMS can be consid- 
ered as another major technique for ultratrace anal- 
ysis. It is — up to now — the only method for 
ultratrace analysis in thin films — a task which is of 
great interest to study the process related contami- 
nation of metallization layers. As an example for 
the application of SIMS, figure 2 shows the ultra- 
trace analysis of the impurity content in molybde- 
num as a function of the number of electron beam 
refining steps. The increase of impurity levels after 
the 6th step is due to the inhomogeneity on the 
cm-scale of the cast molybdenum [11]. 

3. High-Accuracy, Depth Distribution 
Analysis of Dopant Elements in Silicon 

SIMS is presently the most important technique 
for distribution analysis of dopant elements in semi- 
conductors [4,13-17]. For silicon under optimized 
conditions (table 1), detection limits in the range 
between 10' 4 and 3Xl0 13 cm- 3 can be obtained. 
The accuracy of the concentration values in a pro- 
file is on the order of 5% (B) to 20% (P); that of 
the depth scale ~5% (using actual crater depth 
measurements) (see also [13]). The accuracy of 
SIMS profiles can be determined in two different 
ways: 

a) Comparison of profiles quantified with the in- 
ternal calibration ("fluence") method with 
those obtained using external (homogeneous) 
reference materials. Fluence measurements by 
an independent technique like NAA or RBS 
are valuable to check the accuracy of fluence 
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values obtained from measurement of the im- 
plantation current. 

b) Comparison of profiles obtained with differ- 
ent techniques exhibiting independent system- 
atic errors [16,19]. Such comparative analyses 
are of great significance to study the artefacts 
[18] which occur in SIMS depth profiling. 
For all dopant elements, only electrical mea- 
surements [e.g., spreading resistance (SR)] 
can be applied to annealed specimens. There 
is however a systematic deviation between el- 
emental and electrical profiles at concentra- 
tions above the (electrical) solubility limit and 
at the deeper side of the profile, for which SR 
profiles seem to be steeper. Comparison of 
elemental profiles is even more difficult: only 
for 10 B a "reference technique" — NAA exists 
which exhibits the sufficient accuracy — at 
least for annealed specimens [19]. For non-an- 
nealed materials, differences between SIMS 
and NAA depth profiles for B in silicon have 



been observed [20]. For the dopant elements 
As, P, and Sb, chemical etching has to be 
applied for NAA depth profiling. Moderate 
depth resolution and rather large analytical 
errors are encountered. For Sb, RBS can be 
applied for depth profiling. Although the 
method is considered to be largely free from 
systematic errors, it exhibits only a rather 
poor depth resolution and a high detection 
limit. Thus the depth profiles obtained with 
RBS can only be compared with SIMS pro- 
files in a very limited manner. 

Summarizing this discussion of the assessment of 
the accuracy of SIMS depth profiles (which is pre- 
sented in detail in ref. [16]), it must be stated that 
there is a significant lack of independent methods 
to study SIMS artefacts. The most successful ap- 
proach still seems to be the systematic investigation 
of ion-solid interactions, thus creating a broader 
understanding for physical processes like sputter- 
ing and secondary ion formation [21]. 



Table 1. Optimized analytical conditions and practical detection limits for surface distribution analysis 
of dopant elements in silicon with SIMS [4]. ("TT^pr. = practical mass resolution used during measure- 
ment. 
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Up to this point the discussion has been concen- 
trating on "pure" silicon matrices. In industrial 
practice the combined system Si0 2 /Si is usually of 
great interest. This means depth profiling through 
an insulator/semiconductor system with a chang- 
ing matrix. Compensation of charging and matrix 
effects is necessary to achieve quantitative results. 
This problem is especially serious if an element 
must be measured at high mass resolution, like P 
(M/AAf ~4500). New techniques are necessary to 
obtain stable and accurate profiles through the 
layer system SiCVSi [22]: They apply computer- 
ized peak centering routines in switching between 
references mass and analytical ion to eliminate 
magnetic drift and hysteresis effects. Secondly, a 
complete charge compensation is achieved by mea- 
surement of the intensity of the reference ion 



( 30 Si + ) at the steep flank of the energy of distribu- 
tion. This intensity value is extremely sensitive to 
charging: a surface charge of only 1 volt influ- 
ences the reference intensity by 10%. Changes in 
this reference ion intensity are used to monitor the 
process of charge compensation by adjustment of 
the sample potential. The third step — compensa- 
tion of the chemical matrix effect — is performed by 
oxygen saturation of the sample surface during 
analysis. To control the degree of surface satura- 
tion the secondary ion ratio Si0 2 + /Si 2 + , which is 
extremely sensitive towards oxygen coverage of 
the surface, is measured. With this procedure it is 
possible to control the compensation of the chemi- 
cal matrix effect eight times more accurately than 
by monitoring only atomic ions (e.g., the reference 
mass). 



512 



Volume 93, Number 3, May-June 1988 

Journal of Research of the National Bureau of Standards 



Accuracy in Trace Analysis 



These extensive computer based measurement 
techniques finally allow the measurement of the 
distribution of P across thick oxide films into the 
silicon matrix with high accuracy (estimated to be 
~30% for a point on the profile) for determining 
the technically very important assessment of the 
segregation of this dopant element during oxidation 
of silicon in device manufacturing (fig. 3). 

Distribution analysis of dopant elements in semi- 
conductors by SIMS is by far not fully developed 
yet. Due to the increasing demands posed by semi- 
conductor technology new challenges arise, e.g., 
for simultaneous multielement depth profiling, for 
which new computerized measurement approaches 
have to be established [23], or in high lateral resolu- 
tion dopant analysis in devices with finely focused 
ion beams [24], or for high depth resolution profiles 
of very shallow or multilayer structures [25], or 
finally the problem of quantitation through layers 
of changing matrices [26]. 



4. Quantitative Interface Analysis 

Quantitative analysis at interfaces is subject to 
pronounced artefacts and erroneous results due to 
cascade and recoil mixing, changes in sputter rates 
and ionization yields and ion beam induced diffu- 
sion (e.g., for alkali elements) [27], Extensive inves- 
tigations are usually necessary to achieve a 
quantitative result, as shown here in an exemplary 
manner for the study of the distribution of the do- 
pant element Cr between GaAs and the covering 
Si 3 N 4 layer [28]. Figure 4 shows the SIMS profiles 
of i2 Cr + , 28 Si 2+ and 75 As + obtained with oxygen 
primary ions in a GaAs material before and after 
annealing. Cr peaks are observed in GaAs adjacent 
to the interface (A) in both cases and in the SijN 4 
layer (B) after annealing. Ion implantation of Cr 
into a Si 3 N4/GaAs structure also showed these two 
peaks, indicating that peak A was an artefact due to 
chemical yield enhancement by oxygen at the in- 
terface analysis of the annealed sample. Analysis 
with Ar + primary ions showed only peak B, thus 
raising the question about the origin of chemical 
signal enhancement. Modelling of the primary ion 
implantation and sputter process yielded the origin 
of peak B. Due to the higher stopping power of 
GaAs compared to Si 3 N 4 for the primary oxygen 
ions, these are enriched in the near interface layer 
of GaAs when sputtering through the interface 
(enrichment factor ~2). Secondly the sputter rate 
of GaAs is a factor of 3.3 higher than that of Si 3 N 4 . 



Both effects combined lead to a signal increase for 
Cr in the near-interface zone of GaAs of about a 
factor of 10. 

Once peak A was identified as an artefact, the 
real distribution of Cr could be established by fit- 
ting a Gaussian profile to peak B. The accuracy of 
this method was tested by comparison of the 
amount of outdiffused Cr from GaAs with that cor- 
responding to peak B. An agreement for both Cr 
values of 30% was found. 



5. Further Challenges 

Quantitative SIMS has by far not reached its lim- 
its. Due to the increasing necessity for quantitative 
distribution analysis of trace elements particularly 
in high technology materials [29] a strong motiva- 
tion for the development of instrumentation and 
analytical procedures is provided. Some of the 
most significant areas of current and future activi- 
ties are compiled in the following (subjective) se- 
lection: 

i) Analysis of laterally heterogeneous systems, 
e.g., metals with precipitates: reduction of se- 
lective sputter effects by chemical surface re- 
actions to obtain depth profiles with a high 
dynamic range (fig. 5) [30], combination of ion 
microscopy and mass spectrometry to estab- 
lish calibration curves [31]. 

ii) "Monolayer" analysis — quantitative analysis is 
limited by the low number of atoms available 
for a data point, e.g., the detection limit for Al 
in a thin Si0 2 layer on silicon is ca. 
1 ug/g for a consumption of 0. 1 atomic layer 
(e? A =250 \im)~ corresponding to a removal of 
10,000 Al atoms [32]. Laser resonance ioniza- 
tion SIMS employing high transmission TOFs 
promises to extend detection limits to the ng/g 
range or some tens of atoms, respectively [33]. 

iii) Three-dimensional quantitative distribution 
analysis, employing elaborate image process- 
ing techniques. The resistive anode encoder 
based systems [34] seem to be particularly 
promising. Availability of the full 3-D infor- 
mation will allow correlation of all informa- 
tion which can be gained with SIMS— a 
tremendous advantage for (usually complex) 
technical materials. 
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iv) Phase identification by quantitative evaluation 
of cluster ion intensities in the secondary ion 
mass spectrum. Pattern recognition techniques 
have been shown to have great value for se- 
lection of the appropriate features most typi- 
cal for a particular compound or phase [35]. 
Identification of phases (precipitated in a ma- 
trix or at interfaces) is in principle possible 
even if their size is below the spatial resolution 
of SIMS. 

v) Application of chemometric techniques to ex- 
tract and verify information obtained from 
SIMS spectra and profiles, e.g., factor analysis 
[36] or fuzzy theory [37] to evaluate depth 
profiles. 

vi) Compilation of sputter coefficients, ionization 
probabilities, useful yields [38], and relative 
sensitivity factors. These data would be of 
great value for choosing optimal measurement 
parameters and could enable assignment of at 
least semiquantitative values to secondary ion 
intensities for systems where a suitably 
matched reference material is not available. 
Recent investigations on semiconductor [39] 
and refractory metal [10] matrices indicate 
that RFSs can be transferred from one matrix 
to another by the use of scaling factors. 

vii) Development of (further) models for sputter- 
ing and ionization which reflect the physical 
processes occurring under reactive ion beam 
bombardment. 
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Figure 1, Scheme for quantitative multielement ultratrace analysis of refractory metals 
(example: tungsten) with SIMS [9]. 
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Figure 2. Quantitative ultratrace analysis of impurities in molyb- 
denum with SIMS: Relation between impurity concentration 
and numbers of electron beam refining steps [11]. 
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PROCESS PARAMETER: 30 HH Sl0 2 , IMPLANTATION: P— Si, 
50 KEV, 5 , 10 14 CM" 2 , ANNEALING: 1. 1000*C, 10 H1N., 
2 * 3 t HCL,-*-60 NH Si0 2 J 2. lQOD'C, 10 HIN. Nj, 
3. 100 nm CVD Sl0 2 , 900'C, 260 min., N 2 
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Figure 3. Quantitative phosphorus profile measured at high mass resolution {M/ 
AJW =4500) in the SiO^Si system for determination of the segregation coefficient. 
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Figure 4. Depth profiles of 70 nm Si3N,i/GaAs doped with S^xlO^cm - ' Cr. 
Left: before annealing. Right: after annealing. Primary ions: oxygen [28]. 
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Figure S. Quantitative depth profiles of C, O, Al and Si in FeAl 
(0.24), Si(0.25) alloy after annealing in CO-atmosphere [30]. 
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Introduction 

Quantitative electron microprobe analysis has 
been restricted traditionally to application at single 
points or along scan vectors generated by beam or 
stage scans, However, from the very beginning of 
this field, the direct visualization of elemental spa- 
tial distributions by means of x-ray area scans or 
"dot maps" has been a powerful adjunct to quanti- 
tative analysis [1]. A dot map is created by marking 
the beam positions at which characteristic x-rays 
are detected by writing a white dot on a cathode 
ray tube scanned in synchronism with the beam po- 
sition on the specimen. The spatial distribution of a 
selected elemental constituent can be depicted with 
a dot map at the micrometer level of resolution. 
Because of the great inherent value of information 
presented in a visual fashion, dot maps have be- 
come a major operating mode of the electron mi- 
croprobe. However, dot maps are qualitative in 
nature and suffer from several key deficiencies. 
Quantitative information is lost because count rates 
are not recorded; the image information is 
recorded on film and is inflexible for subsequent 
processing; the ability to detect mass concentra- 
tions below 10 weight percent is severely limited 
by time constraints; and the method suffers from 
poor sensitivity to small concentrations changes 



measured against a high average concentration 
level. 



Compositional Mapping 

The development of true compositional mapping 
with the electron microprobe has recently been 
achieved in the microprobe laboratories at the 
National Bureau of Standards working in conjunc- 
tion with the National Institutes of Health [2]. The 
basic methodology consists of recording under 
computer control scan matrix arrays of actual x-ray 
count rates, derived from both wavelength-disper- 
sive and energy-dispersive x-ray spectrometers. 
These x-ray intensity matrices are then subjected, 
on a pixel-by-pixel basis, to complete conventional 
quantitative analysis procedures. These procedures 
include all of the usual steps, such as deadtime cor- 
rection, background correction, standardization 
against intensities measured on pure elements or 
simple compounds, and matrix factor correction to 
account for effects due to electron scattering, x-ray 
absorption, and secondary fluorescence [3]. In ad- 
dition to these conventional corrections, specific 
corrections must be applied to scanned images to 
compensate for the effect of defocussing of the 
wavelength-dispersive spectrometers which arises 
when the beam is scanned off the optic axis of the 
instrument. This defocus correction can be applied 
either mechanically by scanning the specimen stage 
or by rocking the diffraction crystal in synchro- 
nism with the scan. Alternatively, the correction 
can be generated mathematically by making use of 
the geometric equivalence between spatial scan- 
ning on the specimen with a fixed spectrometer dif- 
fraction condition and scanning the x-ray peak 
with the spectrometer while the beam is fixed on 
the specimen. This equivalence can be used to cal- 
culate the spatial dependence of the defocussing at 
any magnification from a spectrometer scan of the 
peak [4]. 

The resulting matrices of concentration values 
can be displayed as images in which the gray level 
or color encoding is based not merely upon the 
x-ray intensities but rather upon the actual concen- 
trations of the elemental constituents. An example 
of a quantitative compositional map of a zinc diffu- 
sion zone at the grain boundaries of a polycrys- 
talline copper sample is shown in figure 1 [5]. This 
image, which was collected in a matrix scan of 9 
hours duration, depicts concentrations as low as 0.2 
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